hep-ph/0506188 
YITP-05-33 



CP violation in the B — > K£ + £ decay 



T. M. Aliev" *, S. Rai Choudhury fe t, A. S. Cornell^ and Naveen Gaur fe § 

"■Department of Physics, Middle East Technical University , 06531, Ankara, Turkey, 

Institute of Physics, Baku, Azerbaijan. 
b Department of Physics & Astrophysics, University of Delhi, Delhi - 110 007, India. 
c Yukawa Institute for Theoretical Physics, Kyoto University, Kyoto 606-8502, Japan. 



Abstract 

Standard Model (SM) CP asymmetries in B — ► K£ + £~ are expected to be very small. 
This feature could help in the understanding of new physics scenarios which predict 
the existence of CP odd phases in various Wilson coefficients. In this paper we have 
analyzed the B — > K£ + £~ decay in beyond the SM scenarios where the Wilson coeffi- 
cients have new CP odd phases. The sensitivity of the CP asymmetries on these new 
weak phases is discussed. 
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1 Introduction 



One of the key ingredients in the Standard Model (SM) is CP violation, which can be described 
by the Cabibbo-Kobayashi-Maskawa (CKM) matrix[l]. However, even with this description we 
still have an incomplete picture concerning the origin of CP violation in the SM. The exploitation 
of CP violation from the theoretical and experimental sides of physics is very exciting, as it may 
open a window to the existence of new physics beyond the SM. Note that the existence of CP 
violation is a well established fact in K[2] and B[3] meson systems. 

In order to study the sources of CP violation it is promising to consider those observables 
which are sensitive to the possible CP phases. For example, CP asymmetries in decay widths and 
lepton polarization asymmetries, such as explored in references [4-7]. 

One of the promising directions for measuring CP violation is the analysis of rare semi-leptonic 
decays. From the experimental perspective the exclusive decay modes, such as B — > K£ + £~ and 
B — > K*£ + £~, are easy to measure. Two years ago the Belle [8] and BaBar[9] collaborations 
announced the following results for the branching ratios for the B — > K£ + £~ and B — > K*£ + £~ 
decays; 



The analysis for study of possible CP violation in B — > K*£ + £~ was done in earlier works 
[10, 11]. The goal of our present work is to similarly study the possible CP violation asymmetry 
in the exclusive B — > K£ + £~ decay using the most general form of the effective Hamiltonian, 
including all possible forms of interactions. Such an analysis will be useful for comparisons with 
experimental results, as the inclusive modes are generally hard to measure. Note that the CP 
violation in the decay B — > K£ + £~ is induced by the b — > s£ + £~ transition, which in the SM 
is practically equal to zero. This is due to the CKM factors V u bV* s being negligible, with the 
result that the unitarity condition produces only an overall phase factor in the matrix element. 
Therefore the CP asymmetry is strongly suppressed. As such, any deviation from zero for the CP 
asymmetry would be an indication of new physics. 

This paper is organized as follows. In section 2, using the most general form of the effective 
Hamiltonian, we derive the matrix element of the B — > K£ + £~ decay in terms of the B — > K 
transition form-factors. We also derive in this section the general analytic expression for the CP 
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violating asymmetry. Section 3 contains our numerical analysis of the CP violating asymmetries 
together with our conclusions. 



2 The matrix element for the B —> K£ + £ decay 

In this section we calculate the matrix element for the B K£ + £~ decay, which is governed by 
the b — > s£ + £~ transition at the quark level. The most general form of the effective Hamiltonian, 
for the b — > s£ + £~ transition (in terms of the twelve model independent four- Fermi interactions) 
can be written in the following form [12]; 



oiGf - 
V2n 



C SL ( sia^Lb) hn + C BR (sia^Rb 



+cl ot L (sl-tM hin L + cf R {s L iM l Rl n R + c RL {s Rll j> R ) l Ll n L 

+C RR (s R ^b R ) £ R ^£ R + C LRLR (s L b R ) £ L £ R + C RLLR (s R b L ) £ L £ R 
+Clrrl (sl b R ) £ R £ L + C RLRL (s R b L ) £ R £ L + Crsa^Ua^l 
+iC TE e^ al3 sa^bea al3 e 
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where L/R — |(1 =p 75), the CVs are the Wilson coefficients of the four- Fermi interactions and 
% — (Pb — Vk)h = ip+ + P-)n is the momentum transfer. Among the twelve Wilson coefficients 
several already exist in the SM. For example, the first two terms with coefficients C,sl and C R r 
describe the penguin operators, where in the SM these coefficients are equal to —2m s C^ and 
—Im^C^. The next four terms in Eq.(l) are the vector type interactions with coefficients C^l, 
C|°pj, C R l and C RR . Two of these vector interactions, Cj£[ and C^, also exist in the SM with 
the form (Cg^ — C\q) and (Cg^ + C\q). Therefore we can say that the coefficients C^l and C^ R 
describe the sum of the contributions from the SM and the new physics, where they can be written 
as; 

itot _ r< e ff 



ritot /nre/j r 1 _J_ C 

^LL — °9 ~~ °10 + <^LL , 

Citot /~ieff 1 f-i \ H 

U LR — 9 + °10 + ^LR ■ 



The terms with coefficients Clrlr, Crllr, Clrrl and C R l R l describe the scalar type interactions. 
The last two terms, with the coefficients Ct and Cte, describe the tensor type interactions. 

Now that we have the effective Hamiltonian, describing the b — > s£ + £~ decay at a scale /i ~ m B , 
we can write down the matrix elements for the B — > K£ + £" decay. The matrix element for this 
decay can be obtained by sandwiching the effective Hamiltonian between B and K meson states; 
which are parameterized in terms of form-factors which depend on the momentum transfer squared, 
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Q 2 = {Pb — Pk) 2 = (p+ — P-) 2 - It follows from Eq.(l) that in order to calculate the amplitude of 
the B — > K£ + £~ decay the following matrix elements are required; 

(K \s llA b\ B) , (K \si<r^<fb\ B) , (K \sb\ B) , (K \sa^b\ B) . 

These matrix elements are defined as follows [14-16]; 

(K(p K )\s^b\B(p B )) = f + 



(Pb + Pk) 



ml 



m 



K 



m. 



m 



K 



(2) 



(K(p K ) \sa^b\ B(p B )) = -i ■ (p B + Pk)^u - q^PB + Pk)v ■ (3) 

frig + m K L J 

Note that the finiteness of Eq.(l) at q 2 = is guaranteed by assuming that /+(0) = /o(0). 

The matrix elements {K(pk) \sia^ v q u b\ B(pb)) and (K \sb\ B) can be obtained from Eqs.(2) 
and (3) by multiplying both sides of these equations by g M and using the equations of motion, we 
get; 

(%) \sb\ B( PB )} = f m| ~^ , (4) 

m b - m s 



{K(p K ) \sia^ v q v b\ B(p B )) = f (p B + p K )^q 2 - q^(m 2 B - m\) . (5) 

m B + rax i 



As we have already mentioned the form-factors entering Eqs.(2)-(5) represent the hadronization 
process, where in order to calculate these form-factors information about the nonperturbative 
region of QCD is required. Therefore for the estimation of the form-factors to be reliable a 
nonperturbative approach is needed. Among the nonperturbative approaches the QCD sum rule 
[13] is more predictive in studying the properties of hadrons. The form-factors appearing in the 
B — > K transition are computed in the framework of the three point QCD sum rules [14] and in 
the light cone QCD sum rules[15, 16]. We will use the result of the work in [16] where radiative 
corrections to the leading twist wave functions and SU(3) breaking effects are taken into account. 
As a result the form-factors are parameterized in the following way [16]; 

/j(g2) = 1 - q 2 jm\ + (l-qi/mj) 2 ' (6) 



r 2 



where 1 = + or T, and 
with mi = 5.41GeV and the other parameters as given in Table 1. 



fo(q 2 ) = , i^r, (7) 

1 - q 2 /m 2 ftt 
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n 


T2 


™% 


u 


0.162 


0.173 




fo 


0. 


0.33 


37.46 


h 


0.161 


0.198 





Table 1: The parameters for the form-factors of the B —> K transition as given in [16]. 



Using the definition of the form factors given in Eqs.(2)-(5) we arrive at the following matrix 
element for the B — > K£ + £~ decay; 



M(B^K£ + r) = ^-V tb V ts 



+ ££Q + h 5 £N 



A(Pb+Pk)h + Bq^ 

+hV [C{PB + Pk)„ + Dq^ 
+4£a" u £ (-iG) (p B + Pk)^ - (p B + Pk)v% 

+A£a a(3 £ H (p B + p K )^v - (pb + Pk)v% 



The functions entering Eq.(8) are defined as; 



A 

B 

C 
D 



(C L °l + C L ° R + C RL + C RR )f + + 2{C B r + C S l) 
{Cll + C l °r + Crl + C RR )f- — 2(C B r + C S l) 

K^LR + <^RR - <^LL - <^RL)J+, 
\ U LR + U RR ~ U LL ~~ U RL)J-, 



h 



It 

(m B + m K )q 2 



(m, 



Q — fo — —{ClRLR + C R LLR + Clrrl + CrLRl), 

m b - m s 



N = fo- 
G 



m n — m 



m b - m s 



— (Clrlr + Crllr — Clrrl — C R l R l), 



H 



™b + rriR 

Cte 
mB + m K 



fr, 



(8) 



(9) 



where 



/- = (fo ~ U 



m%-m' K 
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From Eq. (8) it follows that the difference from the SM is due to the last four terms only, namely 
the scalar and tensor type interactions. For an analysis of the CP asymmetry it is necessary to 
compute the differential decay width for B — > Ki + i~ . From the expression of the matrix element 
given in Eq.(8) we calculate the following result for the dilepton invariant mass spectrum; 

f g (B - KIH-) = \V tb V*\ 2 AV'(1, r K , s)vA(s) , (10) 

where A(l,f^,s) = 1 + f\ + s 2 — 2tk — 2s — 2rxs, s = q 2 /m 2 B , fx = rn 2 K /m 2 B , = m^/mB, 
v = a/1 — 4m| / s is the final lepton velocity, and A(s) is; 

Atyi 2 r 

A = ^Re - 96Xm B m e (AG*) + 2Am 2 B rh]{\ - r K )(CD*) + \2m B m t {\ - r K )(CN*) 

3 

+l2m 2 B m 2 e s \D\ 2 + 3s \N\ 2 + 12m B m £ s(DA^*) + 256Am^sw 2 \H\ 2 + Am|(3 - v 2 ) \A\ 2 
+s3sv 2 \Q\ 2 + 64Am^s(3 - 2v 2 ) \G\ 2 + m|{2A - (1 - v 2 )[2\ - 3(1 - r K ) 2 ]} \C\ 2 ] . (11) 

As we have already mentioned, our goal in this work is the study of possible CP violating 
asymmetries beyond the SM in the B — > K£ + £~~ decay; at this point we shall briefly remind the 
reader of the situation in the SM . In the SM the Cg Wilson coefficient is the only one to have strong 
and weak phases. Strong phases arise from the short distance effects and resonances whereas the 
weak phase comes from the CKM elements. The remaining two coefficients, CV and Cio, are 
strictly real within the SM. From the parameterization of the form-factors it follows that they 
are inherently real and thus the imaginary parts in the functions in Eq.(ll) can come only from 
the Wilson coefficients in Eq.(l). By strong and weak phases we mean the phases which are CP 
even and odd respectively. In other words we shall consider the picture where CP violating effects 
due to the short distance dynamics are parameterized by the Wilson coefficients. In principle all 
Wilson coefficients can have nonzero strong and weak phases . 

In general the amplitude for B — > K has the general form [11]; 

A(B -> K) = e^Aie** 1 + e i<t>2 A 2 e i&2 , (12) 

where the strong phases are labeled as <5's and the weak phases by 0's. As noted above the strong 
phases are CP even, whereas weak phases are odd under CP. Thus we arrive at an amplitude for 
the conjugated process, B — > K, from Eq.(12); 

A(B -> K) = e-^Mie** 1 + e"^ 2 A 2 e i5 \ (13) 

where the amplitudes of the decay rate of particle and anti-particle can be defined by the CP 
asymmetry (in the decay rate) as; 

\A\ 2 -\A\ 2 = -2A 1 A 2 Sin(<l> 1 -(l> 2 )Sin(5 1 -5 2 ) 
CP \A\ 2 + \A\ 2 A 2 + 2A 1 A 2 Cos(<t) 1 -(f) 2 )Cos(5 l -5 2 ) + A 2 2 [ ' 
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Note that from the above expression we observe that in order to have CP asymmetry we should 
have both strong and weak phases in the amplitude; where the strong phases are provided by 
Cg . In the SM the weak phases for the b — > si + i~ transition are negligible and hence the CP 
asymmetry for processes based on the quark level transitions, b — > s£ + £~, are highly suppressed. 
We will now consider the CP asymmetry in the decay width which is defined as; 

^(B - K£ + £-) - ^(B - Ktt-) 
Acp(q 2 ) = § % • (15) 

as as 

Note that one can also have a CP asymmetry from the Forward-Backward (FB) asymmetry [23]. 
However, in our present case the FB asymmetry for B — > Ki + i~ vanishes within the SM. 

We shall now consider the minimal extension of these Wilson coefficients. In this approach 
we shall assume that the Wilson coefficients corresponding to scalar and tensor type interactions 
vanish identically (of course in the general case we can consider all Wilson coefficients with an 
arbitrary weak phase). For scalar type operators which emerge in Supersymmetric (SUSY) models 
and two Higgs Doublet Models (2HDM) this assumption is justified when we have electrons or 
muons in the final state. The reason being that in SUSY and 2HDM these operators originate 
from an Higgs exchange which results in Wilson coefficients which are proportional to me, and 
hence negligible for i — e, 

The Wilson coefficients for the dipole operator obeys; 

C BR =-2C e 7 ff m b , C SL = -2C e 7 ff m s , (16) 

with 

C 7 e// = |CV e// |exp(20 7 ), 

where <p 7 is an arbitrary phase and it is not constrained by the already observed branching ratio 
Br(B -> K*j). 

Regarding the appearance of the new weak phase in Cio we feel that a few words are in order. 
One of the possible discrepancies between the experimental results[17] and the theoretical predic- 
tion for B — > TtK (from the B — > rnv data) can be resolved, as proposed in [18], by introducing a 
complex phase in the Wilson coefficient C w = Cf M exp(i0i O ). In this prescription the weak phase 
given to C 10 does not effect the CP asymmetry in B — > Kl + l~ . 

We will assume that the Wilson coefficients C RL and C RR also have weak phases, that is; 

C RL = \C RL \ exp(i(fr RL ), 

C RR = \C RR \exp(i(f) RR ). (17) 
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The Wilson coefficient Cg^(m b ,q 2 ) has a finite phase, where, in order to better appreciate 
this, we write its explicit phase content as; 

C e 9 ff (m b ) = C 9 (m 6 )(l + (§) ) + Y SD (m b , s) + Y LD (m b , s) , (18) 



7T 

where Cg(m b ) = 4.334. Here ou (s) represents the 0(a s ) corrections coming from the four quark 



operator Oq [19]; 



2 2 4 _ _ N 2, _ , /a _ 5 + 4s 



- (-) = ~ 3^ (5) - 3 In (5) In (1 - 5) - 3ir ^ iy In (1 - 5) 

2g(1 + g)(1 - 2g) ln( g)+ , 5 + 9 *" 6 * ■ (19) 



3(l-s) 2 (l + 2s) W 3(l-s)(l + 2s) " 

In Eq.(18) >so and Y LD represent, respectively, the short and long distance contributions to 
the four quark operators CW,...,6 [19, 20]. Here Y SD can be obtained by a perturbative calculation; 

Y SD (m b , s)=g (m c , s) [3d + C 2 + 3C 3 + C 4 + 3C 5 + C 6 ] 

-^(M) [4C3 + 4C4 + 3C5 + C6] 

(0, s) [C 3 + 3C 4 ] + [3C 3 + C 4 + 3C 5 + C 6 ] 

"W^[3Ci + C7 2 ][ j7 (0,S)- ff (7h c ,S)] , (20) 

Vts V tb 

where the loop function g(m q ,s) represents the loops of quarks with mass m q at the dilepton 
invariant mass s. This function develops absorptive parts for dilepton energies s = 4m 2 ; 



g (m q , s) = --\nm q + — + -y q - - (2 + y q ) - y q \ 



x 



(e(l-y ff ) ( ln 1 + V^^ _ CT j + e( j/g -l)2arctan ^= |, (21) 



where m ? = m q /m b and y g = 4m^/s. Therefore, due to the extension of the absorptive parts of 
g (m q , s) we see that the strong phases come from Y SD . In particular one notices that the terms 
proportional to g (0, s) have a non-vanishing imaginary part independent of the dilepton invariant 
mass. 

In addition to these perturbative contributions the cc loops can excite low-lying charmonium 
states *p(ls), ■ ■ ■ : tp(6s) whose contributions are represented by Y LD [21]; 

Yld (m b , s) = -J-^ C<°> - §^ [3C 3 + C 4 + 3C 5 + C,\ } 
« I VuVfb V ts V tb J 

,V r , (Ml -Sml-iMy.Ty) ' { ' 
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where Ki is a phenomenological parameter taken here to be 2.3 so as to produce the correct 
branching ratio of Br(B -> J/ipK* -> = Br (B -> J/ipK*)Br(J/ip -> [5], and = 

3Ci + C2 + 3C3 + C4 + 3C5 + C*6 = 0.362. Contrary to Y50 the long-distance contribution in Yk£) 
has both weak and strong phases. The weak phases follow from the CKM elements whereas the 
strong phases come from the s values for which the i-th charmonium states are on shell. Therefore, 
the Wilson coefficient (m b ) has both weak and strong phases already in the SM. 

In this sense the Wilson coefficients Cy (ra&) and Cio(mb) can not develop any strong phase, 
and thus, 07 and 0io should necessarily originate from physics beyond the SM. As such the phases 
of 07 and 0io can be chosen to have a purely weak character. 

3 Numerical analysis 

In this section we present our numerical results for the asymmetries A C p for the B — > Kfi + \i~ 
decay. Note that the parameters for the hadronic form-factors are taken from Table I. For values of 
the Wilson coefficients in the SM we have used C 3 = 0.011, C A = -0.026, C 5 = 0.007, C 6 = -0.031, 
C^* = —0.313, C 9 = 4.344, and C w = —4.664. For further numerical analysis the values of the 
new Wilson coefficients are needed, where we have varied them in the range — |Cio I < Cx < |Cio|- 
The experimental value of the branching ratio of the B — > K(K*)i + i~ decays[8, 9] and the bound 
on Br (B — > /i + //~)[22] suggest that this is the right order of magnitude. It should be noted that 
the experimental results lead to strong restrictions on some of the Wilson coefficients, namely 
—2 < C LL and C RL < 2.3, while the remaining coefficients vary in the range — j C10 1 < C x < |Cio|- 
For the remaining parameters we take m b = 4.8GeV, m c = 1.35GeV, m B = 5.28GeV and m K = 
0.496GeV. 

For the kinematical interval the dilepton invariant mass is 4mj < q 2 < (mg — m^) 2 where 
the J /ip family of resonances can be excited. The dominant contribution comes from the three 
low-lying resonances J/ip : ip' in the interval 8GeV 2 < q 2 < 14.5GeV 2 . In order to minimize the 
hadronic uncertainties we will discard this subinterval in the analysis below by dividing the q 2 
region in to low and high dilepton mass intervals; 

Region I : 4m 2 e < q 2 < 8 GeV 2 , 

Region II : 14.5 GeV 2 < q 2 < (m B - m K ) 2 , (23) 

where the contribution of the higher resonances do still exist in the second region. 

As mentioned previously, we have analyzed the case where there are four weak phases: 7 , 0i O , 
4>rl and 0^. 

In fig. 1 we have presented A C p in the 07-g 2 plane for the B — > Kfx + ^~ decay for Region I and 
Region II respectively. In Region I the CP asymmetry is practically independent of q 2 , becoming 
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Figure 1: Plot of the CP asymmetry in B — > Kfi + fi~ as a function of the phase 0/6*7(07) an d the 
dilepton invariant mass. The left plot is for Region - I and the right plot is for Region - II. The 
other Wilsons are taken to have their SM values. 

maximal in the value for CP violation at <p 7 = n/2. In Region II, however, the q 2 dependence is 
comparatively enhanced as the dominance of the dipole coefficient is now reduced. Aside from 
this our figures suggest that the CP asymmetry in Region II is four times larger than in Region 
I, and this confirms our earlier expectation. 

Since the CP asymmetry is dependent on q 2 and the new weak phases there can appear 
some difficulties. The dependence of one of the variables, for example q 2 , can be removed by 
integrating over q 2 in the allowed practical kinematical region, where the averaged asymmetries 
could be measured more easily experimentally. Therefore we shall now discuss only averaged CP 
asymmetries, which we define in the following way. That is, our averaging procedure is defined 
by; 




where Ri means Region I or II. 

We now depict in fig. 2 the 7 dependence of the averaged asymmetries (A C p). From this 
figure it can be observed that the average (CP) asymmetry can attain values of 3%. Differences 
from zero of any value of (A C p) would be an unambiguous indication of the existence physics 
beyond the SM. 

In fig. 3 we have plotted the dependence of CP asymmetry on the dilepton invariant mass 
and <p RL . In fig. 5 we have shown the same kind of plot but for 4>rr. We have also shown the 
correlation of averaged CP asymmetry and the integrated branching ratios. In fig. 4 the variation 
of (Acp) with integrated branching ratio for B — > Kji + ji~ for Crl is shown. In this figure we have 
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Figure 2: The averaged CP asymmetry ({Acp) ) in the B — > K^i + ^ decay, where C-j has a phase. 
The remaining Wilsons are taken to have their SM values. 

used three different values of Crl and have varied the phase ((f) rl) in the range < 4>rl < 2?t. 
All the other Wilsons are taken to have their SM values. In a similar graph, given in fig. 6, we 
have varied Crr. 

In the present work we have studied the sensitivity of the CP violating asymmetry on the new 
weak phases appearing in the Wilson coefficients. We have also observed that the CP asymmetry 
in Region II is 4-5 times larger than that observed in Region I when we consider a weak phase fa. 
This can be understood in that in Region II contributions coming from other operators become 
comparable with the dipole operator Oj, where this operator is dominant in Region I. Having 
obtained the averaged (Acp) asymmetry we obtained a maximal value of approximately 3%. 
Note that an additional weak phase in C\q will not give rise to any CP asymmetry, however, if 
non-standard 5 electroweak operators are considered then the CP asymmetry in the region of high 
dilepton invariant mass can reach a value of up to 10%. 

As stated earlier, we can also, in principle, have weak phases in scalar and pseudo-scalar 
operators. The presence of weak phases in these operators can also substantially effect the CP 
asymmetry. The popular extensions of the SM, such as SUSY and 2HDM, all predict the existence 
of such operators. However, the magnitude of these Wilson coefficients is predicted to be small 
when the lepton £ = e or \i. In the presence of these operators one also gets a non-zero value for 
the FB asymmetry in B — > K£ + £~ . The FB asymmetry could provide another measure of CP 
asymmetry [23] which has not been considered in this work. 

The observation of CP asymmetry in B — > K£ + £~ would not only tell us about the nature of 
weak phases but would also give us an insight in to the structure of the effective Hamiltonian. 
5 by non-standard we mean operators which are not present within the SM 



Region - 1 
Region - II 
Total 
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Figure 3: A plot of the CP asymmetry in B — > Kfi + fi~ as a function of the phase 4>rl o,nd the 
dilepton invariant mass. In this plot we have taken \Crl\ = 2 and with the other Wilsons as 
having their SM values. 

Therefore the measurement of the CP violating asymmetry would provide us with an useful insight 
into the mechanism of CP violation, which in turn would serve as a good test for physics beyond 
the SM. 
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